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We report idealized mode-coupling theory results for the glass transition of ensembles of model 
fullerenes interacting via phenomenological two-body potentials. Transition lines are found for 
Ceo, C70 and C96 in the temperature-density plane. We argue that the observed glass-transition 
behavior is indicative of kinetic arrest that is strongly driven by the inter-particle attraction in 
addition to excluded-volume repulsion. In this respect, these systems differ from most standard 
glass-forming liquids. They feature arrest that occurs at lower densities and that is stronger than 
would be expected for repulsion-dominated hard-sphere-like or Lennard- Jones-like systems. The 
influence of attraction increases with increasing the number of carbon atoms per molecule. However, 
unrealistically large fullerenes would be needed to yield behavior reminiscent of recently investigated 
model colloids with strong short-ranged attraction (glass-glass transitions and logarithmic decay of 
time-correlation functions) . 
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I. INTRODUCTION 

Fullerenes, hollow cages of carbon atoms (Cat), were 
discovered and first synthesized in the mid-ninetcen- 
cighties 0, Q The behaviour of Ceo has been partic- 
ularly thoroughly studied. There are many reasons for 
the interest in these molecules, notably their stability 
and near-spherical structure as single molecules. 

Liquids composed of fullerene molecules have long been 
proposed, but there are two main fundamental issues con- 
cerned with their existence. Molecular dynamics sim- 
ulations (Refs. [E El IE IE 0) predict that individual 
fullerene molecules should remain intact up to temper- 
atures of several thousand K, and arc discharge experi- 
ments indicate that they are stable on at least the 0.1 ms 
timescale. However, when fullerenes are heated in bulk 
to relatively low temperatures (around 1200 K), a dis- 
sociation into forms of amorphous carbon sets in (e.g. 
Refs. H 0, EE EH)- This instability could be due to 
impurities [rj ]. or be a collective effect involving inter- 
fullerene interactions 

Secondly, even if the stability of the high-temperature 
bulk fullerene system is assumed, there has been con- 
siderable disagreement on the characteristics of its equi- 
librium phase diagram, in particular the existence or ab- 
sence of a thermodynamically stable liquid phase (l3lll4j . 
Ceo appears to be a borderline case 0, 0, and the 
gas-liquid spinodal eventually shifts into the metastable 
region below the fluid-crystal transition with increas- 
ing number of carbon atoms per molecule, N 0, Il8j . 
For this reason, model fullerenes have been mentioned 
as analogues of short-range attractive systems 0, |2(| . 
Such systems can be realized in colloidal suspensions 
where added free polymer induces an effective interaction 
whose width and strength can be fine-tuned. They are 
frequently modeled by a square- well system (SWS), i.e., 
a hard-sphere (HS) potential supplemented by a square- 
well attraction of finite range. As this range is decreased, 



the gas-liquid critical temperature drops, until it eventu- 
ally falls into the fluid-crystal coexistence region. Thus 
the metastability of the liquid phase can be taken as a 
signature of 'short-ranged' attractions (see Refs. [2"lH22T| ). 
Indeed, common models for fullerene interaction assume 
LJ-like interactions between the individual carbon atoms, 
which are of short range compared to the large diam- 
eter of the Cn molecule itself. Furthermore, at high 
enough temperatures, the internal degrees of freedom 
of the fullerene molecules should not play a major role, 
since each molecule is able to rotate freely j^E). In this 
way, fullerenes suggest themselves as a molecular system 
that might display the hitherto purely colloidal class of 
phenomena associated with short-ranged interparticle at- 
traction. 

If one is concerned with the kinetic phenomena of the 
(perhaps metastable) fluid, this analogy raises an in- 
teresting question. In experiments on colloid-polymer 
mixtures 0> HE and theoretical studies on the SWS 
or similar models pfjl l27j , qualitatively new dynamical 
features are connected to the short range of the attrac- 
tions. In particular, these systems form glasses; that is, 
amorphously arrested states, whose nature changes as 
the importance of interparticle attraction over the core 
repulsion changes. At very high temperatures, the glasses 
formed from dense fluids are driven almost exclusively by 
the core repulsion - as particles become trapped in cages 
of their neighbors - and are hence qualitatively similar to 
hard-sphere glasses. Addition of attractive interactions 
that are comparable to or bigger in range than the typi- 
cal interparticle separation in this high-density fluid has 
a small effect on the kinetics, since these add up to a more 
or less flat background [2^]. In particular, the glass tran- 
sition in the Lennard- Jones (LJ) model (after all a good 
representation of many molecular liquids) is, according to 
recent theoretical studies |29j . repulsion- driven over the 
entire temperature range. The localization length of the 
particles in such glasses is comparable to 10% of the core 
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diameter, in agreement with the Lindcmann criterion for 
freezing, and little influenced by the range of the inter- 
particle attraction. However, if sufficiently short-ranged 
and strong attractions are present, a qualitatively differ- 
ent glass is found when the thermal energy becomes com- 
parable to or lower than the attraction depth, as bonding 
provides a second particle-trapping mechanism. Conse- 
quently, the range over which particles are localized in 
such an attraction- driven glass scales with the range of 
this attraction and can be an order of magnitude smaller 
than the one predicted for repulsive glass. Furthermore, 
their degree of arrest is much larger; in particular they 
appear mechanically much stiffer than comparable stan- 
dard glasses. The intriguing question then posed by this 
analogy is: do fullerenes, under certain conditions, form 
attractive glasses, or at least ones that display some char- 
acteristics of the colloidal attractive glasses? 

In this paper, we explore theoretically the above anal- 
ogy on the kinetic level between the liquid of model 
fullerenes and short-range- attractive model colloids. Re- 
cent molecular-dynamics (MP) simulations on model 
fullerene liquids |30l l3ll l32 | have shown that there oc- 
curs a positional glass transition in model fullerene liq- 
uids at high density, connected with a freezing-in of the 
molecules' translational degrees of freedom, and not to 
be confused with the rotational glass transition occur- 
ring in fullerene crystals below room temperature |33j |. 
From the MD data, it was concluded that the positional 
glass transition is essentially a repulsion-dominated glass 
transition reminiscent of hard spheres. However, this 
judgement was based solely on the shape of the tran- 
sition line Q g (T). Our study shall provide a theoretical 
basis for the investigation of such fullerene glasses, by 
applying the mode-coupling theory of the glass transi- 
tion (MCT) [2H to the model potentials for fullerenes 
employed in these and similar MD simulations. MCT 
has been very successful in describing the details of the 
repulsive-to-attractive crossover in the SWS [2(| ^ , and 
hence should be well suited to investigate the analogous 
question here. It will allow us to predict not only the 
qualitative shape of the glass transition line, but also the 
strength of the arrest and the form of the time decay of 
density fluctuations near to the transition. The investi- 
gation of kinetic effects might also be interesting for a 
further understanding of the destabilization of fullerenes 
(as far as kinetic cooperative effects are involved). 

We start by giving a brief introduction to the fullerene 
models used in this study (Sec. Hljl. and a sketch of 
MCT and the calculation of glass-transition lines in 
the density-temperature diagram of various fullerenes 
fSec. IIIl|) . Section llVl presents our main results, after 
which a discussion follows in Sec. 



II. FULLERENE MODELS 



the fullerene bulk system by particles interacting with a 
spherically symmetric pair potential. A popular model 
of this kind is due to Girifalco j3f| . This potential can be 
applied to fullerenes formed from any number TV of car- 
bon atoms; however, the assumption of spherical shape 
and the neglect of intramolecular vibrations become more 
and more questionable the larger the molecules. The 
model assumes Lcnnard- Jones interactions with exper- 
imentally determined parameters between the individ- 
ual carbon atoms, which are then integrated over hollow 
spheres representing the buckyball molecules. This yields 
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where r is the centre-centre separation of two fullerene 
molecules of diameter d, and s = r/d. The constants 
a and (3, tabulated in Ref. [35|], are proportional to TV 2 . 
A suitable radius rf(TV) can be inferred from experiment 
and simulatio n fsR ISrl I37I l3?j| . or from the assumption 
d cc \/N [Sfj Eflj . We will use the Girifalco potential 
in our work for TV = 60, 70, and 96, the most common 
low- TV fullerenes. 

Data from ah initio calculations on Ceo can be better 
fitted by a potential due to Pacheco and Ramalho (PR) 
|4l| . One combines a long-range van der Waals expres- 
sion with a short-ranged Morse-potential term, 

V PR {r) = F{r) x M{r) + [1 - F{r)] x W(r) , (2a) 
W(r) = -C 6 /r 6 - C s /r 8 - C w /r w - C 12 /r 12 , (2b) 
M(r) = M exp[r(l - r/d )} ((r(l - r/d ) - 2) , (2c) 

where F(r) describes the cross-over between the two, 



F(r) = (l + exp[(r-p)/<5])- 



(2d) 



Cg and C% may be determined from density functional 
theory, whilst C10, C12, Mq, dn, r, fi, and S are fitting 
parameters tabulated in Ref. |41|. The PR potential is 
somewhat softer than the Girifalco potential and appears 
superior in a number of predictions, but its application 
to larger fullerene molecules is not straightforward. For 
our discussion of qualitative features, the PR potential 
serves to demonstrate to what extent our results depend 
on the precise form of the potential chosen. 

In comparing results for different-sized fullerenes, we 
need to introduce suitable natural scales of energy and 
length. An obvious choice for the energy scale is T* = 
/«BT/e, where e is the magnitude of the potential at its 
minimum. A convenient choice for the unit of length is 
given by the Barker-Henderson effective diameter |42^ . 
letting r* = r/<T e ff with 



For temperatures high enough so that internal modes 
of the Cat molecules become unimportant, one can model 
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FIG. 1: Girifalco potentials for Ceo, C70, and C96 (solid, 
dash-dotted, and dashed lines), in units of reduced temper- 
ature and the Barker-Henderson effective diameter a c g. The 
Lennard- Jones potential is shown as a dotted line. Inset: com- 
parison of the Girifalco (solid line) and the Pacheco-Ramalho 
(dotted) potentials for Ceo- 



Vrop is the repulsive part of the potential only, V lep = 
V(r) for r < r*o, and zero otherwise, where r$ is the 
point for which V(ro) — 0. As T — ► 00, we recover from 
this cr e ff —* he. the hard-core diameter of the Giri- 
falco potential, and, for hard-sphere systems, cr e g = a. 
Note however, that the approach to the hard-core limit 
is extremely slow: for Cgo, Coff ~ 1.821er at the highest 
temperatures considered (T = 50000 K). For later com- 
parison we also note that a purely repulsive potential that 
incorporates all soft-core effects of the Girifalco poten- 
tial may be constructed following the Weeks-Chandler- 
Andersen (WCA) prescription j43j, 



jV(r) + e r < r_ 



r > r_ 



(4) 



Here, r_ is the position of the potential minimum. 
Abramo et al. l3ll l32| use this potential to define a e g 
according to Eq. @ ■ We refer to this as the WC A-cr e ff in 
the following. If one is interested in discussing changes in 
the effective range of the potential with changing N, this 
definition appears disadvantageous, since it includes part 
of the change of r_/ro in the definition of a e g. The ap- 
proach of this effective diameter to the hard-sphere limit 
a is even slower than that of the Barker-Henderson def- 
inition: for C 60 , cr BH ~ 1.827cr at T = 50000 K. The 
advantage of V^ u t lies in the fact that it includes the en- 
tire section of the potential where d r V(r) > , and a 
comparison of MCT results based on V cu t ( r ) and the full 
V(r) will be used in the following to separate out those 
dynamical effects arising from attractive forces. 

The Girifalco potentials for Cgo, C70, and Cg6, are 
shown in Fig. ^ in units of reduced temperature and the 
Barker-Henderson a c g. They show a weak dependence 
of the minimum position on N, r* ~ 1/ (const. + y/~N): 
C96 might therefore be expected to display more evidence 
of attraction-driven glassy behavior than Cgo- Even for 
N = 60, the minimum is significantly more short-ranged 



than in the LJ potential representative for standard liq- 
uids, shown as a dotted line for comparison. Also shown 
(in the inset) is the PR Cgo potential. The softness of the 
repulsive component of the PR interaction means that 
its Barker-Henderson diameter is smaller than that of 
the Girifalco Cgo system, and the potential is shifted to 
higher 7"/<7bf- 

To apply MCT to any of the above potentials, one re- 
quires the equilibrium static structure factor, S(q), cor- 
responding to V(r). The liquid-state theory to approxi- 
mate S(q) is unconnected to MCT. For the Girifalco po- 
tential, we solve the Ornstcin-Zcrnike (OZ) equation nu- 
merically to determine the direct correlation function c(q) 
(and with it S(q) = 1/(1 — pc(q))), using the hypernetted- 
chain (HNC) and Percus-Yevick (PY) closure approxi- 
mations |44j |. These both have a known problem termed 
thermodynamic inconsistency: two equivalent routes to 
calculate the equation of state from them give different 
results. Much effort has been devoted to improving upon 
this deficiency. However, for an application of MCT, this 
issue is usually unimportant. As we shall explain below, 
the MCT integrals are not sensitive to the q — > region 
of c(g), perhaps with the exception of a narrow region 
directly surrounding the spinodal of the model. Instead, 
they rely on a reasonable description of the q > g pC ak 
region, where <7 pca k is the first sharp diffraction peak. 
Similarly, the precise position of the equilibrium phase- 
transition lines, and even whether the approximation pre- 
dicts a stable or only meta-stablc liquid, is not important 
here, other than by typically setting an overall shift in 
the temperature scale. We have chosen HNC and PY 
to indicate boundaries between which many more refined 
closure schemes will vary our results |45| . In particu- 
lar, most of our discussion is based around the PY ap- 
proximation, since we found it to give more reasonable 
results when probing the hard-core part of the poten- 
tial at very high temperatures. The numerical algorithm 
we employ to calculate S(q) is due to Labik et al jifil |: 
we used a wave-vector grid of 4096 points with a cut- 
off Q = 367.6 nm -1 (R = 35 nm for the radial distri- 
bution function g(r)). Our calculated S(q) (using both 
PY and HNC closures) for T = 1900 K, p = 0.74 11m- 3 
compares well (outside the q — > region discussed above) 
with molecular dynamics data for the corresponding state 
point published by Alcmany et al in Ref. |47| . 



III. MODE-COUPLING THEORY OF THE 
GLASS TRANSITION 

Let us present a brief outline of the (idealized) MCT 
description of structural glass transitions. For in-depth 
reviews, the reader is referred to Refs. 0jEEE3- 

MCT aims to describe the dynamical arrest that hap- 
pens in dense liquids: even in the absence of a thermody- 
namic phase transition, the viscosity and, more generally, 
the relaxation times for many time-dependent fluctua- 
tions in the system increase sharply. As the glass tran- 
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sition is reached, the fluctuations are no longer able to 
decay to zero (as it is the case in a liquid), but a certain 
amount remains "frozen" in the still amorphous system. 
The theory describes this as a feedback mechanism driven 
by slow density fluctuations. 

The most convincing evidence for the MCT scenario 
has come from colloidal systems, although the theory 
deals with both colloidal and molecular liquids alike. In 
the latter, however, the MCT transition is "avoided" due 
to additional relaxation paths called hopping processes 
and not captured in the theory. The MCT transition line 
T c (p), or cquivalently p c (T), can still be inferred from ex- 
perimental data by means of scaling laws the theory pro- 
vides for the dynamics close to its idealized transition. 
However, due to the presence of hopping, the calomet- 
ric glass transition T g will be lower, so that T g < T c or 
equivalently p g > p c . The MCT transition point still in- 
dicates a change in transport mechanism from liquid-like 
to solid-like, although not a complete arrest. In this pa- 
per, we are concerned not with a precise determination 
of the MCT transition point, but rather with its quali- 
tative changes upon varying control parameters such as 
temperature T or density p, and with the spatial struc- 
ture of the resulting amorphous solid. These features of 
the glass are usually captured quite well by MCT 

The structure of the glass shall be characterized by 
its Dcbye- Waller factor (also called form factor or non- 
crgodicity parameter, NEP), f c (q): it quantifies the de- 
gree of arrest of density fluctuations at a wave vector 
q in the glass. It is given by f c (q) = lim t _ > . 00 (f> c (q, t), 
the long-time limit of the normalized density autocorre- 
lation function, 4>(q,t) = (g(q, t)*g(q, 0))/(|f>(g, 0)| 2 ), at 

the transition. Here, g(q,t) = J2k=i ex P[*9 -r fcM] are the 
(Fourier-transformed) density fluctuations, r^.(t) are the 
individual particle positions, and (•) denotes a canonical 
average. (\g(q,0)\ 2 ) = S(q) is the static structure fac- 
tor of the system. Even on the liquid side of the glass 
transition and in presence of hopping processes, f c (q) 
can be determined from the time-dependent density cor- 
relation function. <fi{q,t) (as a function of logt) close 
to T c develops a two-step process consisting of a "fast" 
decay towards its "plateau value" f c (q), followed by a 
slow final decay to zero. The time window over which 
4>(q, t) stays close to its plateau extends as T c is reached, 
and similar plateaus occur in many dynamical quanti- 
ties. The mean-squared displacement (MSD) of a tracer 
particle for example exhibits a plateau that is directly 
connected to the localization length r s of this tracer. 
Also the width of the f c (q)-vcTSus-q curve gives an in- 
dication of the inverse localization length in the glass, 
since it closely follows the plateau value for the tagged- 
particle density correlator, f s ' c (q) ~ exp[— q 2 r 2 }, where 
the latter approximation holds for small q. Attractive 
glasses are characterized by a significantly reduced local- 
ization length, and consequently exhibit f c (q)-\eis\iB-q 
curves that follow a much broader envelope than that 
found in standard hard-sphere or LJ glasses. 



Within MCT, the equation determining f c {q) reads 



/(g) 
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to be evaluated at the transition point, m is the MCT 
memory kernel, 

m[f] = £ I £^V(q, k)f(k)f(p) , (6a) 



2q 2 J (2tt) s 

with p = q — k. The V(q, k) constitute the coupling con- 
stants of the theory, and they are determined completely 
in terms of the equilibrium static structure of the system; 
more precisely in terms of the direct and the triplet cor- 
relation function. Lacking a reasonable expression for 
the latter, one usually drops it (approximating three- 
point averages by a convolution approximation; this is 
not connected to the neglect of three-body terms in the 
potential). One then gets 



V(q,k) = S(q)S(k)S(p) 



k 



c(k) 



q p 



c(p) 



(6b) 



In this sense, S(q) is the sole input to MCT. In partic- 
ular, temperature effects enter only through their effect 
on S(q). The fc-integration in Eq. (|6a|) is responsible for 
suppressing finite variations in the form of c{q — > 0), i.e., 
the influence of thermodynamic effects. 

Eq. (JSJ in general has many solutions f(q); the NEP 
f(q) is the largest positive, real of these. As the con- 
trol parameters and thus the coupling constants V are 
varied smoothly, bifurcations in Eq. JSJ can occur, lead- 
ing to a (generally) non-smooth jump for /(g) = to 
f(q) > 0. These bifurcations are the MCT glass transi- 
tions. They are found by a bisection search in p (at any 
given fixed T), repeatedly solving Eq. JSJ numerically un- 
til the transition point p c (T) is located with the desired 
accuracy. To this end, wave vectors are discretized to 
a grid of M points with spacing Aq. The results pre- 
sented below are typically obtained for M = 315 and 
Aq = 0.4, with values as large as M = 600 (Aq = 0.04) 
for the SWS. S(q) was obtained as explained above, on 
a grid with spacing Ap = 0.01 nm~ 3 for each tempera- 
ture and interpolated linearly between those points; since 
S(q) varies smoothly, this interpolation leads to negligi- 
ble errors. The results were checked for representative 
points with finer q (Aq = 0.2) and p (Ap = 0.001 nm~ 3 ) 
grids, and with larger ranges (M = 420). 

Each point on the MCT transition line is character- 
ized by a real number 1/2 < A < 1, called the ex- 
ponent parameter, that can be calculated knowing the 
coupling constants and f c (q)- For the standard glass 
transition there holds A < 1, and the decay of the 
density correlation function in the liquid is, asymptot- 
ically close to the transition, governed by two power 
laws: <f>(q, t) — f c (q) ~ t~ a for the relaxation towards the 
plateau, and 4>(q,t) — f c (q) ~ — t b for the initial decay 
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from this plateau. The exponents a and b are determined 
by A, 



p [nm ] p [nm ] 

0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 
2500 F I 



r(i 



T(l - 2a) 



r(i + b) 2 
r(i + 26) 



(7) 



The case A — * 1 (corresponding to a — ► and b — ► 0) 
signals the breakdown of this power-law asymptotic ex- 
pansion for cj)(q,t), and points with A = 1 are termed 
higher-order glass transitions. In their vicinity, the 
asymptotic shape of correlation functions is better de- 
scribed by powers of logf |5(j. Such higher-order sin- 
gularities have been found in the short-range attrac- 
tive colloidal systems discussed above. There, they in- 
dicate a discontinuous cross-over between two qualita- 
tively distinct types of glass - the one driven by repulsive 
caging, and the one driven by attractive bonding. The 
higher-order transition point vanishes for larger attrac- 
tion ranges, 6 > 5 C , with S c ps 0.045 in the SWS. But 
even for 8 larger than but close to 5 C , there remains a re- 
gion along the glass-transition line T c {p) where A attains 
a maximum close to unity and where precursors of the 
higher-order transition have been found experimentally 
or in simulations [24l I2BI IBH 15^ . Hence the approach 
of A to unity upon lowering the temperature in a system 
with attractive interactions can be taken as indicative for 
the crossover to an attraction-driven glass. 



IV. RESULTS 

Let us start by discussing the MCT results for the (ide- 
alized) glass-transition lines in the temperature-density 
plane, T c {p). Fig. [2] shows the transition lines for Cgo, 
C70, and C96. Results are shown using both the HNC 
(filled squares) and the PY (filled circles) closures to the 
OZ equation, in order to indicate the degree of uncer- 
tainty imposed on the MCT results by different approx- 
imations to S(q). The respective no-solution boundaries 
for the two closures are also shown (HNC: open squares; 
PY: open circles); they can be taken as a rough approxi- 
mation to the gas-liquid spinodal. Note that the binodal, 
as well as the freezing and melting lines, being thermo- 
dynamic rather than kinetic in origin, are not shown. 

There is a small but noticeable difference between the 
glass transition curves calculated within the two differ- 
ent closure approximations. Since HNC and PY can be 
seen as limiting cases for a number of more refined clo- 
sure schemes |45j , we take the difference between the two 
transition lines to be an indication of that error contri- 
bution to T c (p) that is not inherent to MCT. Note that 
even if the difference is small in terms of the critical den- 
sity p c , it can still lead to a significant uncertainty in 
determining T c , as the lines run nearly vertically in the 
T-versus-p plots shown. The PY-MCT transition is con- 
sistently shifted towards lower coupling strength (smaller 
p and higher T) compared to the HNC-MCT line. This 
is different from the estimated spinodal regions, where 
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FIG. 2: MCT-glass transition lines (filled symbols and lines 
in panel d) for C60, C70, and C96 models: in panels (a) to 
(c), solid (dashed) lines with circle (square) symbols are the 
results for the Girifalco potential using the PY (HNC) clo- 
sure for S(q). Diamond symbols with dotted lines correspond 
to the PR potential for Ceo- Lines with open symbols indi- 
cate the loci of points below which no solution for S(q) exists 
within the OZ closure (indication of the gas-liquid spinodal). 
Crosses are MD simulation results from Ref. 1301 . Panel (d) 
shows the PY and MD (Ref. I3(il0 results in terms of dimen- 
sionless quantities p* and T* defined according to the BH 
(left) and WCA (right) prescriptions (see text for details). 
The thin (thick) arrows in panel (d) indicate the T — > 00 
transition point for hard spheres as calculated from MCT (de- 
termined from experiment), {6/n)p* — 0.516 (0.58) 



PY generally estimates lower T than HNC. This under- 
lines that different aspects of the equilibrium structure 
are responsible for the two phenomena. 

Fig. |3 also exhibits results for Cgo based on the PR 
potential together with the HNC closure approximation 
[diamonds and dotted lines in panel (a)]. This consider- 
ably softer potential displaces the MCT glass-transition 
line to even higher densities. The uncertainty coming 
from the potential modeling is much greater than the 
one due to different OZ closures. Also included in the 
figure as crosses are some MD simulation results |3tJ for 
p g (T), determined from the nonequilibrium system fol- 
lowing a temperature-quench. As expected (due to the 
cutting off of the MCT glass transition by hopping pro- 
cesses), p g > p c . But the shape of the two lines is roughly 
similar. 

The fullcrcnc diameter and energy scale both change 
upon varying N, leading to a known shift of the phase 
diagram towards lower p and higher T 0, ; the MCT 
transition line follows this trend as expected. Panel (d) 
of Fig.|5]is an attempt to scale out these broad effects: it 
shows the transition lines in reduced units as discussed 
above, using both the BH and WCA definitions of the ef- 
fective diameter a e g. Results for C70 are omitted to avoid 
overcrowding. The difference between the two rescalings 
is marked - use of the WCA-cr e ff shifts the glass transi- 
tion line to higher p* for all T*. This can be understood 
by recalling that the WCA definition includes a section 
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of repulsive potential (V(r) < 0, d r V(r) > 0) omitted 
in the BH approach, leading to a higher er c ff. However, 
the two representations of our results have two major 
features in common. Firstly, both show a clear bending 
across of the glass transition lines to lower p* at lower 
T* , even after scaling out (through a B s) of the growth in 
the effective size of the repulsive core as the temperature 
is decreased. Secondly, the transition lines for C96 are at 
markedly lower p* than those for Cgo at lower T* , even 
though the different fullerene lines are close (and cross if 
BH-er cff is used) at high T* . Both these points are sug- 
gestive of the enhancement of arrest through (relatively) 
short-ranged attraction in these systems. 

However, the form of the glass transition lines at high 
T* suggests that we might encounter problems in using 
these rescaled plots to make firm conclusions about the 
nature of the glassy behavior. As T* — > 00, one expects 
the potential to become more and more HS-like (since the 
Girifalco potential has a hard-core excluded volume con- 
tribution), and hence the transition line to approach the 
HS value, p* c = 0.516(6/tt) H3 for the MCT-PY calcula- 
tion. Similarly, the simulation values for p* are expected 
to approach the experimentally-determined value for the 
HS glass transition, p* as 0.58(6/7r). Both asymptotic 
values arc indicated in panel (d) of Fig. [21 as vertical ar- 
rows. An ideal choice of effective diameter would collapse 
the glass transition lines to the HS p* (making them verti- 
cal) as soon as the temperature became large compares to 
the potential depth. However, both the definitions used 
here fail to produce any collapse of our MCT results over 
the temperature ranges considered here, despite the fact 
that our highest T* correspond to T w 50000 K. In ad- 
dition, both rescalings move the glass transition lines at 
high T* to higher p* than the HS value - a result of the 
slow approach of these cr e ff to the HS diameter a referred 
to in Section Ull 

Note that in Refs. using the WCA-tr c ff, it was 

concluded that p* f» 0.574(6/7r) is essentially equal to the 
hard-sphere transition value and independent of T* , and 
hence that the glass formation was entirely repulsion- 
driven. However, given the lack of a clear-cut physical 
reason for choosing a particular er ff, one needs to be 
careful in drawing conclusions solely from the numeri- 
cal values of p* defined in terms of such a o c s over a 
limited temperature range. In the following, we will try 
to find find signatures of attraction- or repulsion-driven 
glassy behavior that are independent of the choice of ref- 
erence length-scale. Specifically, we will study the dif- 
ference between the MCT results computed with the full 
V(r) and those using the purely repulsive V cu t employed 
in the definition of the WCA-cr e ff. The exponent pa- 
rameter A (Eqn. UJ and the non-ergodicity parameter 
f c {q) = lim^oo 4> c (q, t) (given by Eqn. [HJ will also be 
calculated, predicting, respectively, the form of the re- 
laxation of density fluctuations and the strength of the 
arrest. 

The first of these calculations is aimed at distinguish- 
ing the influence of the attractive part of the Girifalco 
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FIG. 3: Glass transition lines from MCT with the PY approx- 
imation for Ceo (left panel) and C96 (right panel). The solid 
lines reproduce the transition lines from Fig. [5] The dashed 
lines correspond to the analogous results where the Girifalco 
potential has been replaced by its purely repulsive version, 
V C ut(r), Eq. (J1J. Insets: solid lines are the full-potential tran- 
sition lines, and dashed lines are those obtained from cutting 
of S(q) at low q in the MCT integral, see text for details. 



potential on the glass transition from that of the repul- 
sion. The MCT transition lines are computed for the 
cut-off Girifalco potential, Eq. |0J, which excludes any 
attraction but retains the functional shape of the repul- 
sive core. A similar comparison of the LJ model with 
its cut-off variant yields MCT transitions that are al- 
most indistinguishable from each other |29j |. indicating 
that purely attractive effects are negligible for all tem- 
peratures. 

As shown in Fig. [3] for Cqo and C96, the situation is 
different in the Girifalco potential. The glass transition 
obtained from the repulsive core of the potential only 
(dashed lines) only follows the full-potential result at 
very high temperatures. It shows a small deviation in 
the intermediate-temperature region, and below T* f=a 1 
it deviates significantly: it bends over to much lower den- 
sities, indicating that for T* < 1, attraction does matter 
for the glass transition. 

This bend is a genuine kinetic, attraction-driven ef- 
fect, and in particular not associated with the proximity 
of the spinodal. To confirm this, we plot (insets of Fig.|3J 
a comparison of the full-potential transition line with 
one obtained from a model that does have attraction, 
but no spinodal effects entering into MCT: the dashed 
lines in the insets of Fig. were obtained by cutting off 
the S(q) determined from the full Girifalco potential at 
q < <7pcak/3, thereby eliminating the sharp increase in 
S(q — > 0) signalling the approach to the spinodal. The 
resulting MCT transition is close to the one including the 
spinodal, indicating that the latter plays no important 
role, and in particular cannot be responsible for the dis- 
cussed deviation of the attractive system from the purely 
repulsive one. Note that the small shift visible in the in- 
sets of the figure is due to the details of our cutoff proce- 
dure which overestimates the MCT coupling constants. 
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FIG. 4: MCT exponent parameter A evaluated along the 
glass transition line, for Ceo (solid lines) and C96 (dashed), as 
function of reduced temperature T*. The upper set of curves 
corresponds to the full potential, while the lower corresponds 
to the cutoff potential, Eq. Dotted lines are results evalu- 
ated from the full potential without the spinodal low-g region. 
Inset: A(T) plotted as a function of unsealed temperature. 



Based on the above results we can conjecture that 
the Girifalco-model glass transition in the experimentally 
accessible temperature range shows strong attraction- 
induced effects. There exists a region at relatively low 
densities, 0.8 < p* < 0.95, where the Girifalco model 
displays a glassy region that would not be there for a 
purely repulsive mechanism. We note further that, also 
in Fig. |2| a "relative reentry" can be seen - there is a 
region 1 < T* < 10 where the purely repulsive system 
glassifies as lower densities than the one including attrac- 
tion. Again, one can allude to the SWS here, where this 
phenomenon is observed for 5 < 0.1c, albeit the reentry 
there is with respect to the hard-sphere transition value. 
In our case, the reentry is only relative to the repulsive- 
transition line; that the latter is not at constant p* indi- 
cates that the Barker-Henderson effective diameter does 
not account completely for all soft-core effects. 

If there is indeed a cross-over from a repulsive glass (at 
extremely high T) to an attraction-affected one, the MCT 
exponent parameter A should, according to the square- 
well analogy, show a corresponding change as a func- 
tion of T* , indicating the cross-over region by a peak. 
In Fig. 0] we plot the T*-dependence of A for the Ceo 
and C96 potentials. Indeed, A(T*) rises significantly for 
T* < 1 and displays a peak at T* ss 0.5 in both models. 
This peak is, however, quite close to the region where the 
glass transition line terminates at the spinodal, although 
it is in no way connected to it. The last statement is again 
seen from the X(T*) corresponding to the model with a 
\ow-q cutoff in S(q) (see discussion of Fig. [3J being close 
to the full-S'(g) calculated one (the small rise in A being 
due to the overestimation of the MCT coupling constants 
referred to earlier). In contrast, the purely repulsive sys- 
tems according to Eq. |@J display a value A « 0.73 that 
is essentially independent of T* and consistent with the 
values one gets for the Lennard-Jones or hard-sphere sys- 
tems. Again, this clearly demonstrates the influence of 




FIG. 5: Nonergodicity parameters f c (q) at the MCT transi- 
tion as functions of rescaled wave vector, q* — qcr^g. Upper 
panel: Ceo (solid lines) and C96 (dashed); upper curves cor- 
respond to T* = 0.5, lower curves to T 2 * = 2 (PY approxi- 
mation). The dotted lines show the corresponding results for 
the Lennard-Jones system. Lower panel: f c (q) for the square- 
well system, with attraction range 5 = 0.10 (solid lines) and 
S = 0.08 (dashed), at temperatures T* as 0.22 and T 2 * = 0.5 
(MSA approximation). 



attraction at sufficiently low temperatures. A however 
remains bounded by 0.85 (0.86) for Cgo (Cge), hence no 
higher-order glass-transition singularity is predicted for 
these models. The peak in A(T*) increases somewhat 
with increasing N, i.e., with decreasing attraction range, 
so that one can conjecture the existence of a higher-order 
singularity in giant fullerene systems. But unrealistically 
large TV would need to be considered for this. 

As explained in Sec. 1 1 1 1 1 A determines the exponents 
for the asymptotic description of the time-dependent re- 
laxation functions. To the "standard" value A = 0.73 
correspond exponents a = 0.591 and b = 0.314. In the 
Ceo model, the values drop to a w 0.396 and b f» 0.250 
at the maximum in A; for C96 we similarly get a w 0.381 
and b f» 0.244. However, these figures are still signifi- 
cantly different from zero, so that logarithmic decay is 
probably not observable in fullerene glass formers. 

As mentioned above, attractive and repulsive glasses 
can be distinguished by their localization length, and 
similarly by the shape of the f c (q) quantifying the de- 
gree of arrest. Figure [5] shows these plotted as functions 
of rescaled wave vector, q* = qcr^ s , for T* = 0.5 (roughly 
corresponding to T w 1500 K) and T* = 2. A compar- 
ison with the corresponding LJ result (shown as dotted 
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lines) reveals that the Girifalco model predicts a higher 
degree of arrest than that occurring in the Lennard-Jones 
liquid: the f c (q*)-vcisns-q* curves for Gn are systemat- 
ically above the LJ ones. In particular, there are two 
noteworthy trends exhibited by these curves, one upon 
changing temperature, the other upon changing N. 

For the LJ system, f c (q*) changes relatively little with 
T* , except for the low-g region as the spinodal is ap- 
proached. This increase in the q — > value is connected 
with the divergence of compressibility there. For the Gir- 
ifalco models, f c (q*) changes rather more with tempera- 
ture. Again, the fact that attraction plays an important 
role in the formation of fullercnc glasses whereas it does 
not in LJ glasses, manifests itself here. 

Shortening the attraction range with increasing TV, we 
find the f c (q*) to increase at fixed T* . indicating that 
particle caging is enhanced. Extracting a measure of 
the localization length from these curves, one concludes 
that this length shrinks with shrinking width of attrac- 
tion. However, from the half- width of the / c (<7*)-versus- 
q* curve one still estimates localization lengths of the 
order of 5% for the systems studied here, rendering them 
intermediate between the LJ reference system and the 
truly short-ranged SWS. To highlight the connection to 
this SWS, we show in the lower panel of Fig. [5] similar 
f c {q*) for two such systems, S = 0.1 and 5 = 0.08. This 
leads to localization lengths that are approximately those 
observed in the Girifalco system. Also, the values of A are 
similar: A w 0.81 for 5 = 0.08 and A w 0.79 for 5 = 0.10 
at T*, to be compared with A w 0.83 for Ceo (see Fig.0J. 

Again, as T* is lowered, a cross-over occurs from rel- 
atively narrow f c (q*)-vevsus-q* curves to significantly 
wider ones, signalling the crossover from a repulsive to 
an attraction-driven glass. Note that the numerical val- 
ues for T* < T 2 * cannot be compared immediately, since 
the SWS data has been obtained in the mean-spherical 
approximation (MSA) to S(q), which results in an intrin- 
sically different T*-scale [22]]. In the SWS data shown, 
the crossover is more evident than it is in the Girifalco 
model. Note that neither of the two SWS are short- 
ranged enough to exhibit a higher-order glass transition. 

The comparison with the SWS leads quite naturally to 
the question how to define an effective attraction range 
(5 c ff for smooth potentials such as the Girifalco one. One 
such definition has been proposed by Noro and Frenkel 
|54| and is based on an extended law of corresponding 
states: one compares systems with smooth interactions to 
a SWS at a corresponding state (in terms of rescaled den- 
sity and temperature), for which also the second virial co- 
efficient relative to that of a HS system, is matched. 
Since B\ depends on the attraction range, this gives a 
prescription to determine an 5 e s- From the Noro/Frcnkel 
mapping we read off <5 c ff ~ 0.14 for the T* m 0.5 discussed 
in Fig. O However, it appears that the values of A cor- 
respond more closely to a SWS with a narrower range, 
S < 0.10. This can be rationalized easily: the mapping 
by Noro and Frenkel emphasizes a correspondence of £>2 , 
which is argued to take on roughly constant values close 



to the liquid-gas critical point. The physical mechanisms 
responsible for the change in f c (q), on the other hand, are 
quite different from the spinodal effects, as pointed out in 
connection with Figs. |2 and [3J There is no reason to ex- 
pect that a (5 e ff suitable for mapping the spinodal region 
will be the most useful one to map the glass-transition 
region to that of a SWS. 



V. CONCLUSIONS 

We have calculated idealized glass transition lines for a 
number of fullercnc systems with various sizes, using the 
mode-coupling theory of the glass transition and an ef- 
fective pair-potential description. Our results are readily 
testable using molecular-dynamics computer simulation, 
such as performed recently [3fl l3ll 132^ . 

Glass transition lines are found for densities of the or- 
der of 1 nm -3 for Cgo, bending over to slightly lower 
densities at temperatures below the gas-liquid critical 
point. The results for C70 and C96 are qualitatively simi- 
lar, but shifted in the temperature-density plane accord- 
ing to the change in natural energy and length scales 
of the different systems. At least within MCT, these 
glass transitions, being purely kinetic in origin, appear 
completely unrelated to the gas-liquid spinodal or sim- 
ilar thermodynamic questions, and hence to the long- 
standing question whether fullercnc liquids arc stable or 
only metastable with respect to sublimation. Note also 
that a common way to suppress crystallization and hence 
study the metastable liquid in colloidal suspensions is to 
make use of their polydispersity. It might be interest- 
ing whether a similar polydispersity arising in the pro- 
duction of fullcrcncs [55j could play the equivalent role, 
something that has already been indicated in simulation 
studies of binary fullerene mixtures [8^ | . Finally, in the 
observation of our predicted MCT transition line, one 
will need to consider (in analogy to standard molecular 
liquids) the problem of additional (hopping) relaxation 
processes, which may cutoff the MCT transition as such, 
to different extent at different temperatures. However, 
the observability of our predictions should be not worse 
than in other molecular liquids, where MCT has been 
applied with great success at least in a limited region on 
the liquid-side of the transition . 

The glass transitions we have discussed are strongly 
influenced by inter-particle attraction, which is, at least 
within MCT, a clear contrast to standard molecular 
glasses where interactions are well described by Lennard- 
Jones type potentials. This attraction-domination leads 
to an occurrence of the glass transition at lower densities 
than in the purely repulsive system. It manifests itself in 
several ways, in particular through MCT's exponent pa- 
rameter A and the plateau values of the time-dependent 
correlation functions (or the glass form factors). The 
exponent parameter, A > 0.85 at typical temperatures, 
is found to be significantly higher than in LJ glasses 
(A « 0.73), leading to lower exponents a and b for the 
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asymptotic description of the time-dependent relaxation 
in terms of power laws. It should be possible to observe 
this difference, for example by an asymptotic analysis of 
correlation functions measured in MD simulations |4gj . 
Also, the asymptotic form of the divergence of relaxation 
times or the diffusivity close to T c is governed by these 
exponents, r - \T-T C \~7 with 7 = 1/ (2a) +1/(26). Note 
that 7 w 2.44 in the hard-sphere system, while our calcu- 
lations predict 7 ~ 3.46 for the attraction-affected part 
of the fullerene glass transition. The final decay of the 
correlation function is also often fitted with a stretched 
exponential law, <j>{q,t) f=a A(q) exp[— (t/r(q)) 13 (q)], with 
some A(q) < f c (q) and a stretching exponent (3{q) < 1. 
For q — > 00, j3(q) — > b j^j; thus b can be taken as a mea- 
sure of the 'stretching' exhibited by the correlation func- 
tions: low values of b will correspond to a more stretched 
decay. It is conceivable that such an analysis of dynam- 
ical correlation functions will reveal a much broader re- 
laxation spectrum than in usual glasses. Such tests have, 
to our knowledge, not been performed yet. 

As a second signature, the glass form factors f c (q) 
within the Girifalco model are noticeably higher than 
for corresponding LJ states: the fullerene glass is pre- 
dicted to be relatively stiff, featuring relatively high 
plateaus in the dynamical two-step relaxation process 
over a wider wave- vector range than usually observed in 
molecular glasses: in the LJ system, this plateau has ba- 
sically dropped to zero for qu c s sa 30, while we calculate 
f c > 0.1 still at this wave vector for the Girifalco system. 
There, one needs to reach wave vectors qo~ e g ~ 50 before 
the amplitude of the final relaxation process vanishes. 
From the trend observed in f c (q), we argue that for Cn 
glasses, the localization length (a measure of the average 
cage size) should scale with 1/N, i.e., with the inverse 
size of the fullerene molecule and therefore its relative 
attraction range. This scaling is similar to that observed 
in attraction-driven square-well systems, and opposite to 
that in the standard class of repulsive systems. Com- 
parison of our predictions for f c (q) with simulations or 
experiment would require a rescaling of density and 
probably also temperature [5^ as a result of the tendency 
of MCT to underestimate the coupling required to pro- 
duce arrest I5MI53I. I t is possible that the conclusion of 
Abramo et al |30l l3ll| that the glass transition in Ceo is 
hard-sphcrc-like (from the collapse of their data to a ver- 
tical line using the WCA-cr e ff) is a result of the fact that 
the molecular dynamics transition is shifted to lower tem- 



peratures (higher couplings) than that in MCT. Lower 
temperatures and densities (if accessible) might then be 
required to see the effects of attraction. 

Note that f c (q) is usually predicted from MCT with 
much better accuracy than e.g. T c , leading to often quan- 
titatively correct results. Measurements of f c (q) would 
therefore be highly desirable, allowing both a test of 
MCT and an estimate of its error in terms of density 
and temperature. 

Attraction-dominated in the sense discussed above, 
however, does not imply the phenomena of very short- 
ranged strong attraction, as they have been discussed in 
the field of colloid-polymer mixtures. Their hallmarks 
are logarithmic decays in the time-dependent relaxation; 
in this respect, the fullerene systems we have studied 
are only intermediate in terms of their attraction range. 
However, attempting to map the results to those obtained 
for square- well systems can still prove useful, as our com- 
parison of f c {q) results indicates. 

In particular, they can be used to define an effective 
attraction range of the fullerene (or similar) systems. We 
want to emphasize that such a definition is not unique, 
and might well depend on the physical problem one is 
interested in. Previous discussions have argued mainly 
in terms of a nearly vanishing stable-liquid pocket, em- 
phasizing the metastability of the liquid-gas spinodal as 
a signature of short-ranged attractions. If one, however, 
is interested in high-density kinetic phenomena, the un- 
derlying physical mechanisms mediated by the attraction 
have little in common with those in the vicinity of the 
spinodal, and hence a different measure of effective at- 
traction range needs to be used. Such a measure could 
be based on demanding equality between the f c (q), or, 
less stringently, the localization length at the glass tran- 
sition. If one is concerned with glassy dynamics, the 
best measure might be to introduce a "law of correspond- 
ing glasses" based on the MCT exponent parameter, de- 
manding A(T*) = Asws(T 1 *) 8efi{T*)) to determine 5 e s- 
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